Prof.(Dr.)A.S.N.Murty 

Chapter  -6 

Waves  in  the  Atmosphere 

The  mathematical  solution  of  most  problems  of  dynamic  meteorology  is  exceedingly 
difficult  owing  to  the  nonlinear  nature  of  the  terms  in  the  equation  of  motion.  So  many 
of  the  problems  in  meteorology  can  be  solved  only  when  they  are  linearized.  The 
motion  that  is  to  be  studied  can  often  be  treated  as  a small  perturbation  super  imposed 
on  an  undisturbed  state  of  the  atmosphere.  Such  problems  arise,  for  instance,  in  the 
theory  of  the  origin  of  tropical  cyclones.  According  to  V.Bjerknes,  cyclones  develop 
as  small  wave  perturbations  at  the  boundary  between  two  air  masses  of  different 
density  and  velocity.  If  such  a wave  is  unstable,  its  amplitude  increases  and  a cyclone 
develop.  The  undisturbed  motion  here  is  assumed  as  geostrophic  motion. 

6.1.  The  perturbation  equations: 

For  deriving  the  perturbation  equations,  the  following  assumptions  are  made: 

i.  Both  the  total  and  undisturbed  motions  satisfy  the  equation  of  motion. 

ii.  The  perturbation  quantities 


Let  us  consider  the  instantaneous  motion  ( u ) as  the  sum  of  undisturbed  (U)  plus  the 
perturbed  motion  (u’).  Writing  like  this: 

u = U+uf  y = V + v*  w = ft7  + w'  anci  p = P + pl  (gj) 

The  total  equation  of  motion  can  be  written  as 


du  du  ~ du  ~ du  ~ du  dp 

— = 1-  u 1-  v 1-  w — = -a 1-  fv 

dt  dt  dx  dy  dz  dx 

dv  dv  - dv  - dv  ~ dv  dp 

= h U b V h W — — Ct jU 

dt  dt  dx  dy  dz  dy 

dw  dw  ~ dw  ~ dw  ~ dw  dp 

= h u 1- v 1 -w — = —a s 

dt  dt  dx  dy  dz  dz 


(6.2) 


And  the  equation  of  continuity  is 


.(6.3) 


du  dv  dw 

— + — + — = 0 

dx  dy  dz 

According  to  the  assumption  (i)  the  same  equations  (6.2)  and  (6.3)  are  satisfied  for 
undisturbed  motion  also.  Which  means  they  are  obtained  by  simply  replacing  the 
capital  letters. 
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But  to  find  the  perturbed  motion  (6.1)  are  to  be  substituted  in  (6.2)  and  expanded. 

Let  us  consider  first  the  x equation: 

CS  C) 

— (U  + u')  + (U  + u')  — (U  + u')  + (V  + v')  — (U  + u')  + (W  + w')  — (U  + u')  = -a  — (P  + p')  + f(V+v') 
dt  dx  dy  dz  dx 

(6-4) 

For  undisturbed  motion  we  can  write: 


dU  TTdU  wdU  dU  dP  _ 

+ U + V + W = -a  — + fV  (6.5) 

dt  dx  dy  dz  dx 

To  get  perturbed  motion  subtract  (6.5)  from  (6.4)  and  neglect  the  product  of 
perturbations. 

du  du'  du'  du'  ,dU  ,dU  ,dU  1 dp'  ,,  . , . . 

— + U — + V — + W — + u + v + w = — + fv (6.6a) 

dt  dx  dy  dz  dx  dy  dz  p dx 

Similarly  we  can  write  y and  z component  equations  as: 
dv'  dv'  dv'  dr'  dV  dV  dV  1 dp' 

— + U — + V — + W — + u'—  + v'—  + w'—  = — ~ju'  (6.6b) 

dt  dx  dy  dz  dx  dy  dz  p dy 

dw'  T1dw'  TTdw'  TI7  dw'  ,dW  ,dW  ,dW  1 dp'  ^ , 

+ U + V + W + u + v + w = — (6.6c) 

dt  dx  dy  dz  dx  dy  dz  p dz 

The  equation  of  continuity  also  can  be  written  as: 
du'  dv'  dw' 

— + — + = 0 (6.6d) 

dx  dy  dz 

These  equations  are  called  perturbation  equations. 


6.2.  Acoustic  or  sound  waves: 

Sound  waves  are  longitudinal  waves.  The  particle  oscillations  in  longitudinal  waves 
are  parallel  to  the  direction  of  propagation.  Sound  is  propagated  by  the  alternating 
adiabatic  compression  and  expansion  of  the  medium. 

To  obtain  one  dimensional  sound  waves  with  the  help  of  perturbation  method,  let  us 
consider  that  the  waves  propagate  in  a straight  pipe  parallel  to  the  x-axis. 

To  exclude  the  interference  of  transverse  oscillations  ( particle  motions  at  right  angles 
to  the  direction  of  propagation),  let  us  assume  that  meridional  and  vertical  velocities 
are  zero. 

So  v = w = 0 and  u = u(x,t)  and  doesn’t  depend  on  y and  z coordinates. 

Then  the  equations  of  momentum,  continuity  and  thermodynamic  energy  for  adiabatic 
motions  are  as  follows: 


du  1 dp  ^ 

— + — = 0 

dt  p dx 


(6.1) 
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1 dp  + 

p dt 


du  dv  dw 

1 1 

dx  dy  dz 


= 0. 


As  3rd  and  4th  terms  are  zero  for  transverse  oscillations,  equation  (6.2)  reduces  to 

du 


1 dp  + 

p dt 


dx 


= 0 ( 6.2) 


And  the  thermodynamic  energy  equation  can  be  written  as 

— (ln<9)  = 0 (6.3) 

dt 

(Please  note  under  adiabatic  conditions  0 is  conserved).  Equation  (6.3)  can  be 
modified  taking  the  potential  temperature  equation 

P o 


6 = T 


.(6.4) 


Where  po  is  suface  pressure  =1000mb.  We  may  eliminate  0 in  equation  (6.3)  using 
(6.4).  ' 

Taking  logs  in  equation  (6.4) 


ln0  = lnr  + — In 


^1000^ 

l P j 


= lnr-  — \n(p) (6.5) 


Using  P = RpT  equation  (6.5)  can  be  written  as 


f 


In  6 = In  P\ 


V 


R 


'pj 


[ln(i?)  + ln(y>)]  = — In  p - In  p since  1 — — = — = — 

Y CP  cp  r 


Or  differentiating  with  respect  to  ‘t’ 

— ln<9  = — — In  p-—  In p = 0 (6.6) 

dt  y dt  dt 

Eliminating  p between  (6.2)  and  (6.6),  we  can  get 

' du' 


l d , 

In  p 

y dt 


dx 


= 0 


.(6.7) 


Applying  perturbation  technique  using  the  identities 


u(x,t)=u+u',  P(X,t)=P  + P'’  P(x,,)=P+P'  (6.8) 

Substitute  eqn  (6.8)  in  (6.1)  first  after  expansion.  Then  we  can  get 
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( - 

) 1 

r-  i 

f-  \ 

u+  ur 

+ 

u+u' 

dt 

V 

1 

l ) dx1 

V J 

+ 


1 


( 


^ dx 

P+P 


v 


P+P 


= 0, 


.(6.9) 


J 


Substituting  (6.8)  in  (6.7) 
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First  rewrite  eqn  (6.7)  as  — — In  p + 

y dt 


1 1 dp 

— + 

y p dt 


C du^ 
ydxj 


= 0 or  — + yp\ 
dt  1 


ydxj 


^ du^ 
\dxj 

= 0 


= 0 


f du > 


ttf+“i+HiJ=0 ",li" 

Substituting  (6.8)  in  (6.10)  we  can  write 

(p+  p')  + («+ «' ) (p+  p')  + y(p+p')  ^-(u+u')  = 0 (6. 1 1) 

dt  Ox  ox 

For  expanding  eqn  (6.9) , the  Bousinesq  approximation  may  be  adopted: 


1 


1 


p+p 


l + £- 


P 


PJ 


P 


v1 


1 + 


P 


P ) 


P 


P_ 

P ) 


P 


( V^<1) 

p 


The  equation  (6.9)  can  be  written  as 


du  du'  ~ du  du'  ,du  ,du'  l dp  1 dp' 

dt  dt  dx  dx  dx  dx  dx  dx 

P P 

We  know  the  undisturbed  flow  as 


.(6.12) 


.(6.13) 


du  ~ du  1 dp  _ 

— + u — + — = 0 

dt  dx  dx 

P 

Subtracting  (6.13)  from  (6.12)  and  neglecting  the  nonlinear  terms  associated  with 

, . .du  ,du' 

perturbation  terms  like  u — = u — = 0 we  get 

dx  dx 


.(6.14) 


du'  - du'  1 dp' 

— + u — + — = 0 

dt  dx  dx 

P 

Similarly  we  can  write  (6.1 1)  as 
dp'  ~ dp'  ~ du' 

— + U — + YP — = 0 

dt  dx  dx 

To  find  the  solutions  of  (6.14)  and  (6.15)  let  us  assume  there  exists  a wave  of  the 
form 


.(6.15) 


it’  = Ae,k<x~ct>  and  p’ = Bemx~c,)  (6.16) 

where  k = 2nfL  is  the  wave  number.  If  (6.16)  are  the  solutions  of  (6.14)  and  (6.15), 
they  should  be  satisfied  when  they  are  substituted  in  them.  So  differentiate  (6.16) 
and  substitute  in  (6.14)  and  (6.15)  and  simplify. 

Then  the  equation  (6.14)  becomes: 
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f- 


B 


u-c  = 


.(6.17) 


P 


Similarly  (6.15)  becomes: 


B u-c  = -ypA  (6.18) 

v J 

Eliminating  A and  B from  (6.17)  and  (6.18)  we  get 


B 


f-  \ 
u-c 

V J 


= -yp\ 


B 

( - \ 

> u-c 

v J 


f - y 
u-c 

J 


Y_P 

P 


yRT 


Or  C = u±  ^{yRT ) 


( -:P  = RpT ) 


(6.19) 


Equation  (6.19)  says  that  the  speed  of  propagation  relative  to  zonal  current  u is 
■J(yRT) . This  quantity  is  called  the  adiabatic  speed  of  sound.  This  is  also  called 
Lamb  wave. 

The  mean  zonal  velocity  here  plays  only  a role  of  Doppler  shifting  the  sound  wave 
so  that  the  frequency 


tj  = KC  = K{u±  V(^)} 

Corresponding  to  a given  wave  number  K appears  higher  to  an  observer  in  the  down 
stream  from  the  source  than  to  an  observer  in  the  upstream. 

6.3.  Gravity  Waves: 

Consider  a homogeneous  incompressible  fluid.  The  undisturbed  flow  consists  of  a 
constant  horizontal  velocity  U between  a lower,  horizontal  boundary  ( z = 0)  and  an 
upper  free  surface  whose  undisturbed  position  is  z = h.  the  waves  are  assumed  to  have 
infinite  lateral  extent  in  the  y direction  and  thus  are  essentially  in  x-z  direction.  The 
friction  and  coriolis  forces  are  neglected. 

With  the  above  assumptions,  we  can  say  U = const,  V = W = 0 
1 ^ 

And = g such  that  0 <z<h 

p dz 
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Assumptions: 

1.  The  wave  lie  in  the  xz  plane 

2.  The  coriolis  and  frictional  forces  neglected 

3.  The  bottom  of  the  layer  is  a rigid  earth’s  surface 

4.  the  undisturbed  flow  is  zonal  ( U = const,  V=W=0) 

5.  The  waves  have  infinite  lateral  extent 

Due  to  these  assumptions  the  perturbation  equations  turned  out  as 


Then  we  can  write  the  perturbation  equations  (6.6  a,c  and  d)  as 


du'  TT  du' 
— + U — 
at  dx 


1 dp' 

p dx 


(6.20) 


8w'  TTdw' 

+ U 

dt  dx 


\_djS_ 
p dz 


(6.21) 


du'  dw' 

— + = 0 

dx  dz 


(6.22) 


Assume  that  the  solution  of  equations  (6.20)  to  (6.22)  are  as 

u'  = y/(z)eik(x-ct) 

w'  = 0(z)eik(x-ct)  (6.23) 

p'  = S(z)eikix-C,) 


To  solve  the  equations  substitute  (6.23)  in  equations  (6.20)  to  (6.22)  and  on 
simplification  and  apply  the  boundary  conditions: 


a)  The  kinematic  boundary  condition  says  that  the  vertical  velocity  (w)  vanish 
at  the  bottom  (z  ) i.e.  on  z = 0 , w = 0 

b)  The  second  dynamic  boundary  condition  used  is  that  the  total  pressure 


(P+p’)  is  zero  at  the  top  surface.  In  other  words  at  the  free  surface  —(P+  //)  = 0 . 

dt 


Then  we  can  get: 


(U-C)  = ± 


— tanh(kz) 
k 


or  C = U± 


— tanh  (kz) 
k 


.(6.24) 


As  z = h at  the  top  of  the  free  atmosphere,  on  substitution 
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c = u± 


c = u± 


v 


— tanh  (kh) 

k 


, putting  k =2tt/L  we  get 


gL  ( 2nh\ 

— tanh  

2k  ( L J 


.(6.25) 


The  first  term  on  the  right  hand  side  of  equation  (6.25)  is  called  convective 
term  and  the  second  term  is  called  the  dynamic  term.  The  positive  and  negative 
signs  tell  about  the  positive  and  negative  x directions  of  propagation  of  the  wave. 

Two  special  cases  can  be  discussed  here. 

Case  (i) : Deep  water  waves: 

If  h > 0.4  L,  tan  h (2tt1i/L)  1 


So  Equation  (6.25)  becomes  C = U ± 


gL 

2 n 


-,v. 


.(6.26) 


These  waves  are  called  deep  water  or  short  waves. 

Case  (ii):  Shallow  water  waves: 

If  h < L/25,  then  tan  h (2irh/L)  — > 2irh/L 

Then  equation  (6.25)  becomes  C = U ± yfgh  (6.27) 

These  are  also  called  long  waves. 

For  homogeneous  atmosphere,  h = Po/pg  = 7.991  km,  where  Po/p  = RTo.  Then 
equation  (6.27)  turns  out  to  C = U ± y[RT0  . This  is  nothing  but  equation  (6.19) 

which  is  called  Newtonian  speed  of  sound  and  is  equal  to  280  m/s  at  0°C. 

Note  that  the  dynamic  term  of  equation  6.25  is  real  and  so  the  gravity  waves  are 
neutral. 

6.4.  Atmospheric  Rossby  Waves: 


Rossby  waves  are  very  important  for  large  scale  meteorological  processes.  These 
are  called  planetary  waves  because  they  are  very  long  waves  encircling  the  whole  globe. 

In  a barotropic  atmosphere  these  waves  conserve  absolute  vorticity  ('Q  + f)  which 
owes  its  existence  to  the  variation  of  Coriolis  parameter  with  latitude  called  the  p effect. 

Rossby  wave  propagation  can  be  understood  in  a qualitative  manner  by 
considering  a closed  chain  of  fluid  parcels  initially  aligned  along  a circle  of  latitude. 

From  the  potential  vorticity  theorem  [^r(1T7")  = 0],  consider  the  flow  of  fluid  at 

constant  depth  (D)  on  a latitudinal  circle.  On  this  initial  latitudinal  circle  assume  that  £ = 
0 at  initial  time  to.  Let  suppose  at  time  ti,  8y  is  the  meridional  displacement  of  a fluid 
parcel  from  the  initial  latitude.  Then  at  ti  we  can  write  the  absolute  vorticity  as: 

((  + f)  t-i  = (£  + f)ta 
(Od  + (/)ti  = (0*  + (f)to 
(Oti  + (f)ti  = GO  to  since  (Oto  =0 

Or  (Oti  = -(f)  Cl  + (f)co  = -{(/)  t!  -GOto]  = - (^)  Sy  =p5y  

(1) 
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The  minus  sign  shows  that  the  initial  value  is  less  and  final  value  is  more  (as  coriolis 
force  increases  with  increase  of  latitude). 

From  equation  (1),  it  is  evident  that  if  the  chain  of  parcels  is  subject  to  a 
sinusoidal  meridional  displacement  then  the  CAV  is  positive  (cyclonic)  for  a southward 
displacement  and  negative  (anticyclonic)  for  a northward  displacement  as  shown  in 
figure  below.  This  means  the  pattern  of  vorticity  maxima  and  minima  propagates 
westward.  This  westward  propagating  vorticity  field  constitutes  Rossby  waves. 


Fig.  CAV  trajectories  and  westward  propagating  Rossby  waves. 

The  solid  curve  without  arrows  is  the  latitude  circle.  The  dashed  wave 
shows  the  original  perturbation.  The  solid  wave  with  arrow  marks  shows  the  west 
ward  displacement  of  the  pattern  due  to  advection.  The  arrow  marks  on  the  solid 
wave  shows  the  cyclonic  and  anitcyclonic  motions  and  not  the  direction  of  wave 
propagation.  Note  that  the  positive  circle  is  at  a lower  latitude  than  the  negative 
circle  implies  the  wave  is  westward  propagating. 


Derivation  of  Rossby  waves: 

Assumptions: 

1.  The  flow  is  horizontal  (zonal)  and  non  divergent 

2.  The  atmosphere  is  autobarotropic  so  that  basic  mean  current  ( u ) is  constant  with 
height  and  along  meridian. 

3.  The  perturbation  components  u',  v 1 and  p'  are  independent  of  y.  That  is  they  have 
infinite  lateral  extent. 


Consider  the  barotropic  vorticity  equation:  — ((  + f)  = 0 .... 

d f" 

Or  ^ ^ =0 

dr  dr 

On  the  XY-  plane  (z  component  is  not  considered)  we  can  write: 

= 0 


•(2) 


+ v»1f\  + 

\-+ 

&f 

u — 

+ vdM 

.dr 

dx 

.dr 

dx 

(3) 


Since  — — 0 and  u — = 0 , the  equation  (3)  reduc  es  to 

d £r  diX 


8 . 8 8' 


? + [®  |] 


= 0 


(4) 


Applying  the  perturbation  technique,  keeping  the  above  assumptions  in  view: 
u = u + u , v = v'  , i (5) 
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such  that  f ' = 


u = — 


By  ’ 


v'  = ^ 

5.x 


3t/ 

£-y 

where  \|/  is  stream  function. 

s^'=  =(£)-£(-£)  - S+0  =»*♦ 

X’  = >2*>  p> 

Substituting  perturbations  ( eqn  5)  in  equation  (4) 

;i+(a+u.)£  + *±]r +[.'•$.<> 


.(6) 


•(7a) 


c «{'  -Jf'  r flr'  . ffif'  , ffi/ 

Expanding  ^ + « ^ + « ^ v v * 

Substituting  equations  (6)  and  (7)  in  the  above  equation,  we  can  write  as 

w 

Note  the  third  and  fourth  terms  gets  cancelled.  Use  eqn  (7a)  for  ( {')  in  second  term 
and  — = /?.  Then  the  above  equation  (8)  can  be  written  as: 

£(V»  + «|;(V»+  P(£)=0  (9) 


Let  the  solution  of  this  second  order  differential  equation  (9)  is 

i|/  = ^gdks+iy-vt)  (10) 

then  find  and  substitute  in  eqn  (9).  To  find  V2i|/,  double  differentiate  equation 
(10) 

TL=  (ik)  ei('kx+ly-vt}  (10  a) 

Similarly, 

??=  = -^o(fe2)  e^kx+ly~^  since  i2  = -1. 

OX 

. (11) 

Similarly,  = — x|r0  02)  ellfc*+!y  (12) 

Combining  equations  (11)  and  (12),  ^77  + ^77  = — i|/0  edkac+iy-  t-t)  (k2  4-  J2) 

Using  equation  (10),  V2i(/  = - (k2  + 12)  ih  (13) 

Substituting  eqn  (13)  in  (9),  also  use  equation  (10) 

f [-  (k2  + l2)  i|j  ] + u f[-(fe2  + i2)  ih  ]+  ptya(ik)ei(-kx+ly-^  = 0 

C7k  OX 

(14) 

Using  again  equation  (10)  we  can  write  (14)  as 

£[-  (fc2  + l2)  if/Qe‘(k*+iy-rt)  ] + u j-[-  (k2  + l2)  ^aei(kx+ly-  vt)  ] + 

(ik  ) e'Ck.r+iy  -ut)  _ Q 

(15) 

Solving  equation  (15)  by  differentiating  first  term  w.r.t  ‘f  and  second  term  w.r.t  ‘x’  and 
replace  third  term  using  eqn  (10a) 

- (fc2  + i2)  (-  it?)  ih0eiCka+!y“  yt>  + u [-  (fc2  + l2)  ( ik')  g 

= 0 

Again  replace  i|/  using  eqn  (10) 
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( k 2 + l2)  (-  iv)  + u [-  (k2  + l2)  ( ik)  ^ ] + 

— Nli 


«0] 


= 0 

Differentiate  again  w.r.t  ‘x\ 

- (k2  + l2)  (-  iv)  i|y0eiCk*+!>'-vt)  + u [-  (k2  + l2)  ( ik')  ] + 

P (ik)i\i0ei{kx+ly-vd 
= 0 

Divide  throughout  by  ( i)  \|f0el^jr+i3f-1?t^ 

(fc2  + i2)(V)  + u[-(fc2  + i2)(fc)  ]+  /?(fe)  =0 

(fc2  + l2)  ( v — u k ) + /?  fc  = O 


Or  ( v — uk)  — 
P 


fik 


or  (v)=ufe  — /.  r - ^ or  7 = u 


Pk 

P^J 


T77 3T  — c 


.(16) 


°r  B-c  = imn 

In  the  equation  (16),  k is  zonal  wave  number  and  ‘1’  is  the  meridional  wave  number.  If  k 

2fT  1 

= — , where  A is  the  zonal  wave  length  and  l — where  ‘d’  is  the  width  of  the  wave  in 

A.  a 

the  meridian,  we  can  write: 

P P 


u - c = 


(fc2  + I2) 


4tt: 


/a2 


+ i2 


$ Jl  3, 

If  ‘1’  is  zero  (d  = k>),  then  we  can  write  u-  c = —r 

4?7T* 


or  c = u — 


P A 
4 it" 


.(17) 


This  equation  (17)  is  called  Rossby  equation  for  the  speed  of  long  waves  of  infinite 
lateral  extent  (d  = go)  in  the  westerlies  (minus  sign)  over  mean  zonal  flow  (u  ).  In  the 

_ p jp 

equation  (17),  u denotes  the  mean  velocity  and  — r-  is  the  variation  over  the  mean 
flow. 


Thus  the  Rossby  wave  propagates  westward  relative  to  the  mean  zonal  flow  in  the  mid 
tropospheric  height.  Further  the  Rossby  wave  speed  depends  on  the  zonal  (k)  and 
meridional  (1)  wave  numbers.  Therefore  Rossby  waves  are  dispersive  waves  whose  phase 
speeds  increase  with  wave  length. 

For  a typical  midlatitude  synoptic  disturbance  (say  hurricane)  with  zonal  wave 
length  (K)  of  6000  km  and  latitude  width  (d)  of  3000  km,  the  Rossby  wave  speed  (c  ) 
relative  to  the  mean  zonal  flow  (u  ) from  equation  (17)  is  approximately  (-  6 m/s).  This 
means  synoptic  scale  Rossby  waves  move  slowly. 

If  the  frequency  of  Rossby  waves  is  compared  with  other  waves  from  the  table 


given 

Delow: 

S.No 

Type  of  wave 

Waves  per  day 

1. 

Sound  wave 

~W 

2. 

Long  gravitatinal  wave 

104 

3. 

Inertia  & shearing 

lo3 

4. 

cyclone 

101 

10 


5. 


Rossby 


10' 


From  the  table,  in  the  case  of  a cyclone,  the  frequency  is  10  waves  per  day  whereas  in  the 
case  of  Rossby  waves  ( 10  1 = 1/10)  one  wave  occurs  in  10  days. 
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